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Abstract 

Handover rate is one of the most import metrics to instruct mobility management and resource management 
in wireless cellular networks. In the literature, the mathematical expression of handover rate has been derived 
for homogeneous cellular network by both regular hexagon coverage model and stochastic geometry model, but 
there has not been any reliable result for heterogeneous cellular networks (HCNs). Recently, stochastic geometry 
modeling has been shown to model well the real deployment of HCNs and has been extensively used to analyze 
HCNs. In this paper, we give an effective handover analysis for HCNs by stochastic geometry modeling, derive 
the mathematical expression of handover rate by employing an infinitesimal method for a generalized multi¬ 
tier scenario, discuss the result by deriving some meaningful corollaries, and validate the analysis by computer 
simulation with multiple walking models. By our analysis, we find that in HCNs the handover rate is related to 
many factors like the base stations’ densities and transmitting powers, user’s velocity distribution, bias factor, pass 
loss factor and etc. Although our analysis focuses on the scenario of multi-tier HCNs, the analytical framework 
can be easily extended for more complex scenarios, and may shed some light for future study. 
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1 Introduction 

With the dramatically increasing of wireless traffic as well as the population of wireless terminals, the 
traditional homogeneous cellular network cannot provide sufficient bandwidth for all the wireless terminals. 
In response to the capacity challenges, smaller coverage base stations (BSs) are deployed in hotspots to offload 
and have a range of tens of meters to several hundred meters. This brings heterogeneity to traditional cellular 
network and gives birth to the heterogeneous cellular network (HCN). Heterogeneity is expected to be a key 
feature of the next generation of cellular networks, and an essential means for providing higher network 
capacity as well as expanded indoor and cell-edge coverage. In general, HCNs comprise a conventional 
cellular network overlaid with a diverse set of lower-power BSs such as micro cells, femtocells and perhaps 
relay BSs. 

The BSs in different tiers of HCNs (the tiers of BSs are ordered by transmit power) may share the same 
spectra and have different coverage. Handover happens when a user leaves the coverage of its serving BS 
and handover rate is defined as the number of handovers per unit time. The handover in HCNs can be 
divided into two types: horizontal handover and vertical handover. Horizontal handover is the handover 
between two BSs in the same tier and vertical handover is the handover between two BSs in different tiers. 
Compared with horizontal handovers, vertical handovers are more difficult to implement because the HCNs 
may be deployed by different service providers. Thus, to implement vertical handovers, extra communication 
overhand between the the HCNs is essential. Moreover, vertical handovers would lead to additional transfer 
delay, jitter and high risk of dropping, which degrade the service quality. Therefore, handover rate especially 
vertical handover rate is one of the most metrics to instruct the deployment of mobility management and 
resource management. 

The mathematical expression of handover rate in the homogeneous cellular network has been derived by 
the regular hexagon coverage model [?] and stochastic geometry model [?] [?]. But for the heterogeneous 
cellular network, there has not been any reliable and generalized handover expression, due to the randomness 
of the BS positions of HCNs and the different transmitting power of different tiers. In the real deployment, 
the HCN BSs are distributed irregularly, sometimes in an anywhere plug-and-play manner, which results 
in a high level of spatial position randomness. BSs in different tiers have different transmitting power for 
communication, leading to different cell size for different tiers. As a consequence, it is difficult to characterize 
the cell boundaries and to track boundary crossings by UEs (i.e. handovers) in the global networks. Few 
previous works have resolved the above challenges and give a reliable handover rate expression. 

In the recent years, stochastic geometry modeling has shown its admirable ability of modeling the position 
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distribution of HCN BSs [?] [?], and has provided tractable accurate performance bounds for cellular wireless 
networks. Stochastic geometry is a very powerful mathematical and statistical tool for the modeling, analysis 
and design of HCNs with random topologies. For instance, the modeling is employed for the capacity analysis 
of random channel access schemes like ALOHA [?] and carrier sensing multiple access (CSMA) [?], the 
capacity analysis of single and multiple tier cellular networks [?] and the capacity analysis of cognitive-based 
networks [?]. 

Hence, in this paper, we investigate the handover rates including horizontal and vertical handover rates 
by the stochastic geometry modeling. By employing an infinitesimal method, the mathematical expression 
of instantaneous handover rate is derived for a typical moving UE. From the derivation, we find that the 
instantaneous handover rate is related to the instantaneous moving speed, and is independent of the moving 
direction. That means only the moving speed distribution of the memoryless walking model contributes to 
the handover rates and the handover rates can be derived through averaging the instantaneous handover 
rates by the moving speed distribution. Thus, the derived handover rate expressions are applicable for all the 
memoryless walking models. The derived expressions are validated by computer simulation with multiple 
walking models and the impacts of system parameters like BS density, transmit power, moving velocity of 
UE, path loss factor are evaluated. Although our analysis focuses on the scenario of multi-tier HCNs, the 
analytical framework can be easily extended for more complex scenarios. 

2 Downlink System Model 

A fairly general model of HCNs considered in this paper contains N tiers of BSs that are distinguished by 
their spatial densities, transmit powers, path loss exponents and biasing factors. For instance, as shown in 
Figure 1, high-power macrocell BS networks are overlaid with successively denser and lower power picocells 
and femtocells. Macrocell BSs and femtocell BSs can be well modeled by spatial random processes [?] [?]. 
Under this model, the positions of BSs in the nth-tier are modeled according to a homogeneous PPP (Poisson 
point process) <f>„ with intensity A„ in an Euclidean plane. 

Each BS in the nth-tier has the same transmit power P n , and has the same path loss exponents a n > 2, 
n = 1,..., N. Assume that UEs are uniformly distributed in the Euclidean plane with density of f u , and 
the movements of UEs are memoryless and are independent of the distributions of BSs. Memoryless here 
means the current position of a UE is only related to its latest position and is independent of its more 
earlier position, i.e. P[5(to)|<S , (ti), £(£ 2 ), ■ ■ ■ ] = P[S , (to)|<S'(ti)], where P[x] is the probability of x, S(U) is the 
position of a UE at time U, and ti > ti+ 1, i = 0,1,2,- • • . 
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We assume open access which means a user is allowed to access any tier’s BSs. And consider a cell 
association based on maximum biased-received-power (BRP) (termed biased association), where a mobile 
UE is associated with the strongest BS in terms of long-term averaged BRP at the UE. The BRP from the 
jth BS in the nth-tier is P r ,nj that can be given by 

P r ,nj = P n L 0 (Rnj/ro)~ an B n (1) 

where R n j is the distance of the j th BSs in the nth-tier from the origin, Lq is the path loss at the reference 
distance ro (typically about (Air/v)~ 2 for r$ = 1, where v denotes the wavelength). And B n is the bias factor 
of admission [?], that could extends the cell range (or coverage) of the nth-tier by employing B n > 1. The 
considered BRP is a long-term averaged value and fading is averaged out, and so does not include fading. We 
assume that handover happens only when the UE is going across the boundary of the current BS’s coverage, 
which is determined by the long-term averaged BRP and is shown as Figure 1. 


3 Problem Formulation 

Consider a typical UE that is at the origin and is admitted to the fcth BS of the mth-tier initially. As the 
typical UE moves, it may immigrate to other BSs in the same tier or other tiers. So its admission state can 
be depicted as Figure 2, that at time t , the typical UE is admitted to the fcth BS with probability P a ,k(t ), 
or is admitted to other BS with probability P a ,k( 0) — P a ,k(t). Thus, the instantaneous transition rate from 
k state to the k state at time t is 

Hk(t) = ( 2 ) 

Since the movement of the typical UE is a memoryless process and is independent of BSs’ distributions, 
then the instantaneous transition rate is stable and can be given by 

HT = - lim dPa ’ k ® (3) 

k t—s-o dt v ' 


H'"‘- is the instantaneous handover rate of the typical UE indeed, 
with area S can be given by 


N 


N 


Then, the handover rate in a region 


A h = E[£ H^f u S] = ]T (E [H?))fuS (4) 

m =1 m =1 

where E[a:] is the expectation of variable x, and E [H™] is the average handover rate of the typical UE. As the 
handover is assumed to happen at the boundary of BS coverage, not at the boundary of a specified region 
with area S, the average handover rate is independent of the shape of the specified region. In the following 
section, we would derive the arithmetic expression of E[£f™]. 
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4 Derivation of Handover Rate 

In this section, we would derive E [H™] in two steps: firstly, we use an infinitesimal method to derive the 
instantaneous handover rate at time 0 of the typical UE with instantaneous velocity v, then average the 
instantaneous handover rate by the distribution of the velocity v. Thus, the impacts of walking model on 
the handover rate can be reflected by the distribution of the velocity v, i.e. the walking model decides the 
distribution of velocity v, then further affects the average handover rate. So the analysis is applicable for all 
the memoryless walking models with that thought. Note that the moving direction of the typical UE would 
not affect the instantaneous handover rate and we would give the explanation in the following derivation. 

In the following, we would give the position distribution of the associated BS firstly, and then derive the 
handover probability of {P a ,k{ 0) — P a ,k{t)), the instantaneous handover rate ( H ™) and the average handover 
rate (E [H™]) in turn. 


4.1 Position Distribution of Associated BS 

Denote (R n j,9 n j) as the polar coordinate of the jth BS in the nth-tier. Assume that the typical UE is 
admitted to the fcth BS of the mth-tier initially, thus P r , m k > Pr,nj for all n G {1, • • • , IV}. According to 
the max-BRP based association and the BRP definition in equation 0, the distance boundary condition of 
these unassociated BSs can be derived as equation d5j) based on P r ,mk > Pr,nj , 


Pnj ^ ( 


P B 

1 n-^n 

p B 


■)^(Rmk)^ = R 


lb 


( 5 ) 


where R 1 ^ is defined as the distance lower bound of the nth-tier BSs for the clarity of expression. 

According to the distance lower bound of each tier, the probability density function (PDF) of the associ¬ 
ated BS’s distance ( R m k) can be derived by using the null probability of a 2-D Poisson process with density 
A in an area A, which is exp(— \A). By setting A = A„ and A = tt(R 1 ^) 2 for the nth tier, (n=l,- • • ,N) [?], 
we could give the probability density function (PDF) of R m k as 

N 

f(Rmk ) = 2?rA mRmk exp{-7T ^ X n (R (f) 2 } (6) 

n=1 

oo 

where R = Rmk■ Note that / f(R m k)dR m k — "fm is the probability that the typical UE is admitted to a 

o 

mth-tier BS. 

Since BSs are deployed as a PPP (Poisson point process), the 9 m k is uniformly distributed in the range 
of [0, 27 t], and its PDF is given by 

f{9mk) = ^ (7) 
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Since the distributions of 9 mk and R m k are independent of the coordinate axis and the moving direction 


of UE, we can assume that the X axis is along the moving direction of the typical UE at time 0. 


4.2 Derivation of Handover Probability 

When the typical UE moves an infinitesimal distance of r, i.e. the typical UE moves to the point (r, 0) with 
r —► 0, the BRP from the jth BS in the nth-tier is 


jjnew _ 

*r,nj 


RuBtiLq 


( 8 ) 


(V(Rnj COS (0nj)-r) 2 + (R nJ sill(d nj )) 2 )“" 

According to the max-BRP based association, after the typical UE moves to the new point, it is still 
admitted to the primary BS k only when the BRP from BS k is larger than the BRP from anyone else. 
Given the position of BS k. the probability that the typical UE keeps the primary link to the BS k is denoted 
by Pa,k{r\Rmk,O m k) and is given by equation ©. 


N 


N 




Pa,k(r\Rmk,9mk) = II > P%$] = IT P[cos(0 nj -) < ^ 


(9) 


where P[x] denotes the probability of event x, and x n j is defined as 


^ J =r 2 -( ^"p‘ )“» {R 2 mk -2rRmkCOs{e mk )+r 2 )°^ (10) 

And cos (0 n j) < is derived according to P™mk — Pr,nj i n equation 

Hence, according to equation the typical UE would keep its primary link if all the BSs in the nth-tier 

R 2 .-\-x n ■ 

(n=l,- • • ,N) are in the region of cos (9 n j) < ’ or woun b n °t other wise. 

Jl 2 .-\-x n ■ 

We call the region of cos (0 n j) > %rR~ as ^ le bad region of nth-tier BSs, which is shown as Figure 3. 
It means that if there are nth-tier BSs in the nth-tier bad region, the typical UE would immigrate from the 
serving BS to one of those BSs. Thus, P a ,k{r\R k , 9 k ) equals to the probability that no BS is in its bad region 
for all tiers. 

Denote the area of the nth-tier bad region as A rnn (r. R mk , 0 mk ), then we can give the null probability of 
the the PPP in the nth-tier bad region as exp(A n A mra (r, R mk , 0 mk )). Since all the PPPs {4> n }n=i,— ,jv 
are independent, P a , k (r\R k , 6 k ) is the product of those null probabilities and can be given by the following 
equation 

N 

Pa,k(‘r\Rki @k) = exp( A n-A-mn (rjRmkiOmk)) (H) 

n— 1 
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According to cos(-) < 1 and the definition of bad region, the 9 n j boundary conditions of the nth-tier bad 
region can be given as 


1 > cos (9 n j) > 


R-nj + x nj 

2 rRj 


( 12 ) 


Thus, based on the R n j boundary condition in equation ([5]) and 9 n j boundary conditions in equation 
m, we can further derive the boundary conditions of the nth-tier bad region in the Appendix 17.11 and 
give the results as equation (fl3l) . where i} n j is defined as flnj = — 0™", 9™ ax and are the upper 

bound and lower bound of 9 n j , respectively. The shapes of bad regions with different boundary conditions 
in equation m can be depicted as Figure 3. 

As Figure 3 (a) shows, when R< R m k holds, according to the derivation in Aooendix l7.il the range of 


-), arccos( )] f° r a certain R n j. 


2 rR n 
'-Rmk ' 


Rnj is [R 1 *, r + sJT 1 — x n j] and the range of 9 n j is [— arccos(- 2rRn 

As Figure 3 (b) shows, when both > R m k and \9 m k\ < arccos( ~ R IS k ) hold, the ranges of R n j 
and 9 n j are the same as the ranges of R n j and 9 n j in Figure 3 (a). When both > R m k and 
1 9 rn k | > arccos( ) hold, as shown by Figure 3 (c), according to the analysis in Appendix 17.11 the 
range of R n j is [R l n,r + \Jr 2 — x n j\, and the range of 9 n j is [—7r,7 t] if R n j £ [R l ^,—r + \Jr 2 — x n j], or is 
[- arccos( ), arccos()] if Rnj £ [~r + y/r 2 - x nj , r + yV 2 - x nj \. 


' Rnj 6 0, i?„,- = 0, 

. R n j 6 [R l ni r + V r2 - x nj],j>nj = 

Rnj 6 (R l n, r + s/r 2 — Xnj), &nj = 

(13) 

Based on the boundary conditions, the area of the nth-tier bad region ( A mn (r , R m k,9 m k )) can be derived 
as the equation (fl4l) . where 1(c) is the index function that equals 1 if the condition c holds or 0 otherwise. 
The first term in equation (fill) corresponds to the case of ~ R iT k < cos (9 m k) < which is depicted by 
Figure 3 (a) and (b), and the second and third terms correspond to the case of cos (9 m k) < R ih k that is 
depicted by Figure 3 (c). 


2 7T, 

d2 I 

2 arccos 0 •;& -, 


o ! Rnj+ X nj \ 

2 arccos( 2 2 ), 


Rnj 6 [R l n , —r + \Jr 2 — x n j) 


Rnj e [—r + yV 2 — : 


’ + \A 2 - Xnj] 


COS ( 6 mk) > 
cos (6 mk ) < - 


R n 
R mk 


< COS { 6 mk) < 


r+y/r 2 —Xnj 2 

A m „(r,i? mfc ,0 mfc ) = [ / 2 arccos( n ° )RnjdR n j ]F(— ^jnr < cos (9 m k) < t) 

Rib nj n n 

r+y/r 2 -Xnj 2 -r+y/r 2 -x„j 

+ [ / 2 arccos (-^^ ± )RnjdR nj + f 27rR nj dR„j\I(cos(9 m k) < -%ir) 

-r+^/r 2 —Xnj Rl n 

Thus, the handover probability can be obtained as (P a ,k( 0) — ^{R mk ,9mk}Pa,k{i"\R m k, 9 m k))- 


(14) 
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4.3 Handover Rate Derivation 

Since the instantaneous handover rate H™ is the derivative of handover probability, according to equations 
© and ED, the instantaneous handover rate H™ can be derived as equation ED- 


H rn = _ Hm fP^t) = _ l im ^ {Rmk .e mk} [P*Mr\R mk ,8 mk )] 
0 dt at 


t—y 0 

- -HRmk,e mk } [lim- 


t — >o 

=-E { w..> -<•'£-»-> $)] 


dt 


(a) . N N 

df( ‘ 6Xp( ^n-^-mn (^? Rmk ? ^m/c)))] * ^ 

r— n=l n=l 

r'W'/r l:™ dA mn (r,R mk ,e m k) M ... 

_J11 TRmic,<Wl2^ Hm — JJ • V 

n= 1 r—>-u 


(15) 


where n = lim dL is the instantaneous velocity of the typical UE at time t = 0, (a) is obtained according to 

r —>o at 

the expression of P a k(r\Rmk, 0 m k) in equation (fill) , and (b) is obtained due to lim A mn (r, R m k : 9 m k) = 0 

r—>o 
N 

and lim exp(— J2 A„A m7l (r, R mk , 0 mk )) = 1. 

r -*-° n= 1 

Denote H™~ n as the instantaneous handover rate from the fcth BS in the rreth-tier to the BSs in the 


nth-tier. Thus 


N 


hf = J2 h < 


m—n 

k 


(16) 


According to the derivation of H™ and the expression of H™ in equation (fl5l) , we can give H™ 71 as 


Tjm-n in, r\ 1 - dA mn (r, Rmk, 9 m k) , 

H k = E { R mk ,9 mk }[Xn ■ lim - — -] • V 


(17) 


Based on the expression of A mn (r , R m k,9mk) hr equation (ED, lim 


r —>-0 


dAmn (VjRmk ,6 

dr 


is derived in Appendix 


O and the result is given as equation ED- 


lim dAmn(r ’^ mfc,emfc) = I(cos(6» mfe ) < -^^)[-27 t^^ cos(6» mfe )] 

+I(-^rr < cos (9 mk ) < ^ cos (Omk) arccos( 1 |^ cos (9 mk )) + 2 yj(R l £) 2 - (^^ cos(6» m fc)) 2 ] 


(18) 


Theorem 1. The instantaneous handover rate from mth-tier BSs to nth-tier BSs for a UE with instanta¬ 
neous velocity v is and is given by equation m- 


H™~ n = 8A„A m v 


min(l, 

+ 00 R n 

1 1 i, 

0 0 


1-z 2 


'(Ttf) 2 -* 5 


N 




1 -z 2 


}dz{R l n) 2 exp{-7T ) (A* {R l i) 2 ) }dR mk 


i— 1 


(19) 
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Proof. Based on equations (fI71) and (U8[) , HJf 1 n is further derived in Appendix 17.31 □ 

Hence, the average handover arrival rate Ah in equation ([4]) can be given by equation (l20l) . where f v {v) 
is the probability density function of the velocity v in the specified region with area S and is determined by 
the walking model of UEs. 


N N N 

a h = j: (mrus}) = e e mr n f u s] 

m— 1 m —1 n—1 


min(l 

v ’ pit) / 


+oo N N +oo v 

= fvlS f vfy(v)dv E EI 8 AnA m / / [ 

0 m = ln=l 0 0 




/ 1-2:2 1 1 

\J (^) 2 -* 2 V 

1 UP J 


?lb\ 2 C 


N 

■E( 

2=1 


?lb\ 2\ 


( 20 ) 


Similarly, the average handover arrival rate from mth-tier BSs to nth-tier BSs (denoted as A™ n ) can 
be derived as equation (1211) . 


A™~ n = E [H?~ n f u S] 


+oo 


min ( l ,— 777 -) 
+00 R ' b 


= fuS f vf v (v)dv8A n A m f f [ -w 


1-2 2 


(f±rnk \2 

1-^nr) “ 2 


:+■ 


(% 2 _ 2 2 jV (21) 

Wn) 2 ™ p{-^ E {UR^Y)}dR mk 


Although the derived handover rate expressions are not closed-form, these expressions can be efficiently 
computed numerically as opposed to the usual Monte Carlo methods that rely on repeated random sampling 
to compute these results. 


4.4 Discussions of Handover Rates 

In equations of (l20l) and (ED, the general expressions of handover rates have been derived. In this section, 
some corollaries and special cases of handover rates in the stochastic modeling of the HCNs would be given. 

Corollary 1. A™ -ra = A 1 f l ~ rn holds for any m,n € {1, - - - ,IV} ; that is, the forward and reverse handover 
rates between any two tiers are the same. 

Proof. See Appendix 17.41 □ 

Corollary Q] holds under the condition that the UE movements in different tiers are homogeneous and 
UEs are uniformly distributed. Corollary Q] indicates that the mobility between any two tiers would reach 
statical balance over time. 
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Corollary 2. When N = 1, all the BSs are homogeneous and {A n } = A ; the expression of average handover 
rate A h can be further simplified as 

+oo 

X h (N = 1) = -f u S J vf v (v)dv (22) 

o 

Proof. When N = 1, according to equation of (|2H|) . A h can be further derived as equation (l23l) . where (a) is 
obtained due to R = R m k- 


N / a \ +oo +oo 1 

A h(N = 1) = £ (mrfuS}) = fv$ I vfv(v)dv 8A 2 / j2dzRl k exp{-^XRl k )}dR mk 

m—1 0 00 

+oo +oo 

= fuS f vf v (v)dv 16A 2 / R^k exp{-^TXR^jdRmk 
o o 

+oo 

= fuS J vf v (v)dv^Q( 0) 

+oo 

= ^rfuS J vfy{v)dv 

0 


(23) 


□ 


Corollary [2] is consistent with the handover rate expression of homogeneous cellular network given in [?] 


Corollary 3. When {a n } = a, i.e. BSs in all the tiers have the same path loss exponent, the handover rate 
of X™~ n can be simplified as equation j 24 \), 


AJT n (M 

where /3 n = ( )° • 


\ 2A„A m /3 IfuS 

(*)= N ' - 

w( £ Ai/3 4 2 )i-5 

i= 1 



(24) 


Proof. When {a n } = a , according to equation ([5]), R.’f 1 can be simplified as R l £ = ( p"g” )° Rmk = /3 n Rmk- 
Then A™ - " i n equation (1?I1) can be further derived as equation (05l) . 


2A„A m /9 2 /„S 


*(E>i© LS o 

i=l 


' 00 i-— r-rz —— iv 

A r _n ({ a «} = °0 =/« 5 ' / vf v (v)dv8X n X m / / [Jjpzpz+J^-]iz^lR 2 mk ex. p{-7r(£ Ai^ 2 )^}^, 

0 0 0 V i=l 

+oo min(l,/3„) - r— +oo JV 

= fvS f vf v (v)dv8X n X m Pl f / i?^ fc exp{-7r(E Ai/3 2 )i?^ fe }di? mfe 

0 0 V " 0 i=l 

+oo min(l,/9 n ) -_ r— —- 

= fuS f vf v (v)dv8X n XmPn f [J^r]lz -^- 

0 0 V " v 47r(£ Ai/3 2 ) 1 ^ 

i = l 

+oo min(l,/3 n ) - r— - 

/ vf v (v)dv J [J 

n 0 v v 

(25) 
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□ 


Corollary 4. For a UE with constant velocity v, its residence time in a BS coverage of mth-tier (denoted 
by T™) is exponential distributed with average value jfrk, i■ e. the PDF of T™ can be given by 

1 k 

TJm Tjm 

f(T™) = exp(- —T™) (26) 

where 7 m is the probability that a UE is associated with the mth-tier BS and is given by equation m referred 
to [?]. 


y n 

1 im = 27rA m I Bmk GXp{ 7T ^ ( A n [R n ) }dB, 7n k 
i n=l 


(27) 


Proof. When the velocity v is constant, the instantaneous handover rate H™ would keep constant and do not 
change with time. So the transition rate from k state to k in Figure 2 is constant, i.e. — lim —— = 


—is constant for all tn > 0. Thus 
7 m u — 


due to Pa,k(0) = lm- So, 


H rn 

Pa,k(t ) = 7 m exp(- —t) 

r im 


f(T r m ) = 


dt 


= ^exp(-^T r m ) 

7 m 7m r 1 


(28) 

(29) 

□ 


Corollary 5. When the UEs are not uniformly distributed in the whole region and UEs in each tier BS 
coverage are uniformly distributed, the handover rate from mth-tier BS to nth-tier BS in a specified region 
with area S can be given by equation 13(A) . 


K~ n {{Uu} ± u)=mr n fm,uS} 

R 

’ R l J> ■ 


+OO 


+OO 


= fm,uS f vf v (v)dv8X n X m f f l. pr mk 2 

0 0 0 V 


1-z 2 


(-Wjz)2-~2 


■ R l J> 


N 


fes-^(^n) 2 exp{-7T ^2(Xi(R l ^) 2 )}dRmk 

i= 1 


(30) 


where f m ,u is the UE density in the m th tier BS coverage. And the total handover rate in the specified region 

N N 

is A h({fm,u } ^ fu) = J2 J2 A™ _n ({/m,u} 7^ /«)■ 


m— 1 n— 1 


Proof. For a UE admitted to a mth-tier BS, its average handover rate to a nth-tier BS is On the other 


7m 


hand, the average coverage area of mth-tier BSs in the specified region is S"f m . So the total handover rate 


from mth-tier BSs to nth-tier BSs in the specified region is A™ n = E[—r fm,uS^m] = E[.ff™ n fm,uS]. □ 


'fm 
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Corollary 6. When the walking models of UEs in different tier BS coverage are different, the handover rate 
from the mth-tier BS to nth-tier BS in a specified region with area S can be given by equation W- 


a # Mv))=mr n fuS] 


+oo 

/ 

o 




+oo R‘~ b ■ 

I 

0 0 


{ ^) 2_ Z 2 


N 

= fuS J vf m , v (v)dv8\ n \ m f f [J : , +V ' "i-^ —^(^) 2 exp{-7T E(A i{Rf) 2 )}dR mk 


(31) 


where fm,v(v) is the velocity distribution in the m th-tier BS coverage. And the total handover rate in the 

N N 

specified region is \h({fm,v{v)} ± f v (v)) = J2 J2 A ™~ n {{fm,v(v)} ± fv(y)). 


m= 1 n— 1 


Proof. According to the derivation of Hff~ n , HJf~ n is only related to the velocity distribution in the mth- 
tier BS coverage and is independent of the velocity distribution in other tier BS coverage. Thus, Corollary 
[G] is can be directly derived from the derivation of Hff~ n in the former section. □ 


From the above analysis, we could give the general steps for handover rate analysis by stochastic geometry 
modeling: firstly, obtain the PDF of the associated BS’s position, secondly, by the infinitesimal method, get 
the area of the bad region, based on which, then derive the instantaneous handover probability, and derive the 
instantaneous handover rate through taking the derivative, and at last, average the instantaneous handover 
rate by the distribution of UE’s velocity. 


5 Simulation Results 

5.1 Validation of the Analysis 

Now that we have developed the general expression of handover rate for HCNs by stochastic geometry 
modeling, it is important to see how well the analytical results match the computer simulation. Here, we 
consider two types of walking models, i.e. straight-line walking model and RWP walking model. For the 
straight-line walking model, UE would move without changing its moving direction, while for the RWP 
model, UE would change its moving direction for a randomly chosen direction of [0, 2n\ and then keep the 
direction for a randomly chosen duration of [0,100s]. For both models, the velocity is uniformly distributed 
in [ 0 , 2v\. 

In the simulation, we considered a circle region with radius of lOKm, generated BSs of each tier and their 
positions by the PPP with the density of \ n for each try, and generated the UEs with the density of /„. 
For each try, all the UEs move with a specified walking model for 10 4 seconds. We counted the number of 
handover in a specified region with area S = 1 km 2 and gave the handover rates by averaging the results. 
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A total of two tiers are modeled according to PPP in the simulation (N=2). For the simulation, some 
parameters keep constant and these parameters are f u = 100/km 2 , S = 1km 2 , Pi = 1, P 2 = 0.2, Bi = 1, 
cni = 3.5 and Ai = 1/krn 2 . Other parameters would vary in different tries and these parameters are the 
velocity v, the 2nd-tier BS density A 2 and, 2nd-tier BS’s path loss factor 02 and the 2nd-tier BS’s bias factor 

b 2 . 

Figure 4 compares the analytical and experimental handover rates of A/~ 1 , A^ -2 and A 2-2 . In Figure 
4, the dotted lines represent the analytical results, the circles and triangles represent the simulation results 
obtained by straight-line walking model and RWP walking model, respectively, And we set A 2 = 2/krn 2 , 
B 2 = 1 and a 2 = 3.5. In Figure 4, it is observed that the handover rates of A/ -1 , A^ -2 and A 2-2 increase 
linearly with the average velocity of UE. This matches the derived expression of handover rate. In Figure 4, 
we can see that A/” 1 is the largest, A/ -2 takes the second place and A 2-2 is the minimum for any value of 
average velocity. This can be explained as follows. Since Pi = 1 and P 2 = 0.2, the coverage of lst-tier BS is 
much larger than the coverage of 2nd-tier BS. So the lst-tier BS has longer boundary line than the 2nd-tier 
BS. Thus, the total length of the boundary line between two lst-tier BSs is the longest, the total length of 
the boundary line between a lst-tier BS and a 2nd-tier BS is shorter and the total length of boundary line 
between two 2nd-tier BSs is the shortest. Hence, the handover is most likely to happen at the boundary 
between two lst-tier BSs, is less likely to happen at the boundary between a lst-tier BS and a 2nd-tier BS, 
and is least likely to happen at the boundary between two 2nd-tier BSs. From the figure, it can be seen 
that the relative error between the analytical results and the simulation results are less than 3% for both 
walking models. The relative error is most likely brought by the limited simulation time. The good matching 
validates that the analysis is reliable with the variation of UE velocity. In the simulation, we find that the 
straight-line walking model and the RWP walking model without pause time almost have the same handover 
rates when they have the same velocity distribution. 

Figure 5 illustrates the analytical and experimental results of the total handover rate (A h) and the 
handover rate between the two tiers (A^ 2 ) with different 2nd-tier BSs density (A 2 ) for different average UE 
velocity (u). From the figure, we can observe that handover rates increase linearly with the average UE 
velocity, and both A h and A/ -2 increase with the 2nd-tier BS density. This is because larger BS density 
leads to smaller cells and UEs are more likely to move out the smaller cells. The analytical results match 
well with experimental results of straight-line walking model and RWP walking model, with relative error 
less than 3%. This validates that the analysis is reliable with the variation of both BS density and moving 
velocity. 
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Figure 6 demonstrates the analytical and experimental results of handover rates with variations of 2nd- 
tier BS density (A 2 ) and path loss factor ( 02 ). The analytical results match well with the experimental 
results for both walking models with relative error less than 3%. This validates that the analysis is reliable 
with the variations of both BS density and path loss factors. From this figure, we can observe that all the 
handover rates increase with the 2nd-tier BS density. The handover rate of A^ _1 increases with 2nd-tier BS 
path loss factor 02 , Aincreases with firstly, then decreases with a 2 , and A 2 ~ 2 decreases with 02 - This 
is because the coverage of the 2nd-tier BS decreases with an and the coverage of the lst-tier BS increases 
with 02 according to equation ©. Thus, the boundary line of the lst-tier BSs increases with the 02 and the 
boundary line of the 2nd-tier BSs decreases with the 02 - So, the total length of the boundary line between 
two lst-tier BSs increases with 02 , the total length of the boundary line between two 2nd-tier BSs decreases 
with 02 , and the total length of the boundary line between a lst-tier BS and a 2nd-tier BS increases firstly 
with o 2 , then decreases with 02 - It is the reason that leads to the increasing of A^“ , the decreasing of A 2 ~ 2 , 
and the changing of A^ -2 . 

Figure 7 demonstrates the experimental results of the forward and reverse handover rates between the 
two tiers for the RWP walking model with the variations of 2nd-tier BS density, path loss factor and average 
velocity. Figure 7 shows that the forward handover rates (A^ -2 ) match well with reverse handover rates 
(A 2 ” 1 ) between the two tiers with relative error less than 3%. This validates the reliability of Corollary [I] 

Figure 8 demonstrates the analytical and experimental results of residual time distributions with constant 
velocity. The analysis results (denoted by ‘Theo.’ in the figure) is calculated by equation (l26l) . and the 
experimental results (denoted by ‘Simu.’ in the figure) is obtained by the RWP model with constant velocity 
u=5m/s. From the figure, we can see that the experimental results match well with the analysis results, 
which validates the effectiveness of Corollary 01 As the figure illustrated, the residual time in the lst-tier 
BS is much larger than the residual time in the 2nd-tier BS due to the fact that the lst-tier BS coverage is 
much larger than the 2nd-tier BS coverage. As the figure shown, both the residual time in the lst-tier BS 
and in the 2nd-tier BS decrease with the 2nd-tier BS density A 2 , but the decrement of residual time in the 
lst-tier BS is more significant. This can be explained as follows. With the increase of 2nd-tier BS density, 
the average coverage of each tier BS decreases, so the residual time decreases for each tier. On the other 
hand, with the increase of 2nd-tier BS density, the average lst-tier BS coverage decreases more significantly 
because the coverage of each lst-tier BS is much larger than the 2nd-tier BS coverage and the increased 
2nd-tier BSs would be more likely to occupy the lst-tier BSs coverage. Thus, the decrement of residual time 
in the lst-tier BS is more significant. 
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5.2 Effect of Bias Factor 

Figure 9 shows the numerical results of handover rates with the variation of 2nd-tier BS’s bias factor (B 2 ). 
From this figure, we can observe that the total handover rate (A h) decreases with the 2nd-tier BS’s bias 
factor, A^ _1 and A^ -2 decrease with B 2 and A 2 ~ 2 increases with B 2 when B 2 € [1,2], which is a reasonable 
range of B 2 . This can be explained as follows. According to equation EOl) . when R mk equals , A™ -n 
could reach the minimal value, and A™ _n decreases with R mk if R mk < R * l r l ‘ holds. Thus, when B 2 increases 
and Rik < R l 2 holds, R\ k increases and A^ -2 decreases. When B 2 increases, the coverage of lst-tier BS 
decreases and the coverage of 2nd-tier BS increases according to equation m, so the boundary line of lst-tier 
BS decreases and the boundary line of the 2nd-tier BS increases, which lead to the results that A} -1 decreases 
with B 2 and A 2 ~ 2 increases with B 2 . The results demonstrate that we can decrease the total handover rate 
by reasonably adjusting the bias factors. 

6 Conclusion 

In the literature, there has not been any general handover rate derivation for the heterogeneous cellular net¬ 
works. Thus, in this paper, we give a generalized handover analytical framework by employing the stochastic 
geometry modeling for heterogeneous cellular networks, derive the arithmetic expression of handover rate, 
give some meaningful corollories and validate the analysis by computer simulation. The analysis may shed 
some light for future extension and study. 
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7 Appendix 

7.1 Boundary conditions of bad region 

According to the bad region boundary conditions in equation (1121) . we can further derive as follow: 


finj = 0 , 

'dnj = 2 arccos 
'Onj = 27T, 


Rq 7 + x n 

2rR ni 


R n j+ X rij , 

2 rR n j 

1 ^ Rnj+ X n,j ^ 

1 - 2 rR nj - 1 

Rnj+ X ™j ^ _-1 

2 rR nj ' 


(32) 


For the case of — 1 < 


R ni +Xn 


- 2 rR„ 


< 1, we can further derive the inequations of 


Rli+ X n 


2rRni — 


< 1 and 


Rli+Xr, 


2rR n i — 


> -1 
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as equations (l33l) and (1M1) . respectively, 

Rnj - 2 rRnj + Xnj < 0 

= $’ Rnj € [r \J T~ '■i'nj , T T \J T 2 X n j ] 

R^j + 2r R n j + > 0 

=> Rnj £ (-oo,-r— \/r 2 -a;„j)U(^'+\/r 2 -a; n j,-|-oo) 

Hence, we can give the range of as equation (l35l) according to equation J5|) and R n j > 0. 


(33) 


(34) 


min (^n, | - r+ yjr 2 - x nj |) < R nj < r+ ^r 2 - x nj (35) 

As — r + \Jr 2 — x n j > 0 when r —> 0 + , we can further simplify equation (1351) as 

min(i?^,■ -r + ^r 2 -x nj ) < R nj <r+yjr 2 -x nj (36) 

For case of J '" j > 1, similar to equation (l33l) . we can further derive the inequation of ^ I; " J > 1 
as R nJ e (max( J R^-, r + y'r 2 - x nj ), +oo). 

For the case of " J < — 1, similar to equation (13411 . we can further derive the inequation of " 3 < 

-1 as R n j e (R l *j,-r + -^/r 2 - x„j) 

On the other hand, according to the expression of R in equation ([5|) and x n j in equation (fTUl) . R ^ = 
—r + \Jr 2 — x n j holds when r = 0. And the derivative of —r + \]r 2 — x n j at r = 0 is derived as 


lim 

r—>-0 


d(-r+y/r 2 ^c nj ) _ f Bn p n ^Sjh, - cos(9 mk ) 

dr ~\BmPm> n \ sl rnk) " Rmk 


-1 = - 


r‘° 

Rmk 


COS(6» mfe ) - 1 


(37) 


Thus, —r + Jr 1 — x n j increases at r = 0 if cos (Q m k) < or decreases otherwise. So, 

j R l n>-r+^r 2 ~Xnj: COS {Omk)>-^Tr 

I R l n <—r+y/r 2 —x n j, cos (O m k)<-^nr 
Similarly, the following relationships can be derived, 


R l n >r+y/r 2 -x nj , cos(6» mfc )>-^ 
Rn < r + ^T 2 — X n j , < COS (d m k) < 


(38) 


(39) 


Based on the above relationships, the boundary conditions of equation (15121) can be derived as equation 


CE3D- 
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7.2 Derivation of lim 

r—>-0 


d-A-rnn (t *?Rmk i@mk ) 

dr 


For the clarity of expression, we define the expressions of Amliir, Rmk,0mk) and 

A inn 5 Rmk 5 @mk ) in the equation of (l40l) . 


= [ / 2arccos( JV " j ;y" J < cos(6> mfc ) < %f) 

R ib n3 

r+y/r 2 -x n j 

A$n{r,Rmk,0mk) = [ f 2 arccos( ” ) R nj dR nj ] J(cos(6> mfc ) < - ^r) 

-T-\-yJr 2 —x n j 


r+y'r 2 -! 


< R-n.i+ X rt 


Arnnfri Rmki @rrik) — [ J* ^^RrijdR n j\I (cOs($ m /j;) ^ ) 

Based on those expressions and equation da. the following relationship can be obtained, 

A = /l(l) 1 A( 2) , 4(3) 

-nmn ' -^mn ' ^mn 


-r+y/r 2 -i 


(40) 


(41) 


Then lim dAmn(r 'j irnfc ’ 9mfc ) can be derived by deriving lim dAl ‘nAr,R mk ,e mk ) ^ { = 1 2 3 . For simplicity, we 

r—>0 “ r r—>0 ar 

define </j“(r), ip d (r) and f(r,R nj ) as ip u (r) = r + sjr 2 - x nj , (p d (r) = -r + ^r 2 - x nj , and f(r,R nj ) = 
2 R n j arccos( I " 2 1 r ^ ni ), respectively. 


lim ■ 

1 —>0 


^ = 1 ™ a? / 


=hm[/(r,^(r))f^(r)+/ r) f/(r,i?„,)d^]/( 1 ) 


- lim n [ I £f( r , R nj)dRnj]I [1) 

R'n 

ip u (r) 

= lim [ f 


-2 R„ 


R.-\-Xnn 

nj nj \2 


d x . 
dr 


r ^° K b 
R 2 . (^“( 0) 2 
=” 3 lim J^ 1 ) [ f 

r ^° («n ) 2 




-1 


dx 
d r 


1 -- 


l' x nj) 
Ayr^ 

-1 


j ,1 

t 2 2^/y U «J 


(<P U (r)) 2 

= lim f — 1 

r-y0 ( fl !b )2 VAr 4 ~4r 2 x n:i - (y+x nj -2r 2 ) 

= lim - «+^„,- 2 r 2 

r—>0 L r 


lim —/Wf 

_V O 


'nj n 2 

,— —r— 2 r 


( {R\ 


r—> 6 


arccos( vVr4 _ 4r2 ^ . , . _^ 

— -f (1) [~ 2( ~ r ,'”1 cos(d mk ) arccos(^- cos (9 mk )) + 2^/(i?(?) 2 - (^ 


(^-r-2r 2 )-(y+x„ 3 -2r 2 ) rf ^ 

y^4r 4 —4r 2 x n j + ~ ^ ~ 

=j|=|= ■> - yy/4r 4 - 4r 2 x nj - ((i?(?) 2 + - 2r 2 ) 2 ] 

-,1b / TTTThTo 


(iffi) 2 


COS 


i(^mfc )) 2 


where / (1) = < —^- mK/ _ R ,o ,,- v-v ~ o- 

— (i?Jif) 2 , lim jbx n j = 2< ' r 71 ^ cos (0 m k) and the L’Hopital rule. 


(42) 


= /(-%ir < cos(6» mfe ) < %^), and (a) is given by f(r, ip u {r)) = 0, (6) is given due to lim x nj = 

h n r—t 0 

,2 
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(43) 


^ = 1™ I f{r,Rnj)dR nj I {2) 

ip d (r) 


ip u (r) 


ip u (r) 


= lim[/( r ,^( r ))| : ^(r)-/(r,^(r))^ rf (r)+ / j-f R no )dR no ]I^ 

T j / \ 

ip a (r) 

= [27T R l n(-^ COS(6» mfc )+ 1)-2 tt^^- cos(6» mfc )]/ (2) 

where I*- 2 ) = I{cos{O m k ) < — ^Srr), and (c) is obtained by f(r,<p u (r)) = 0, /(r, i/? d (r)) = 27r, lim -d-i^ d {r) = 


r —>■() 


¥>“ W 


Rmk 


equation 


, / did\2 

cos(^ m fc) — 1 and lim J -£pf(r,R n j)dR n j = —2'K y R n) cos (Qmk), which can be derived similarly as 


i —>-0 


tp d (r) 


1 ™ ^ = i™ t I 27T RnjdRnjlW 


V d R) 


i—y 0 


r—m 


Ri b 


= lim \2Tup d {r) dlf ] j( 2 ) 


(44) 


r —>■ 0 


= -2nR%(-^ cos (9 mk ) + 1 )/ (2) 


which is obtained due to lim ip d (r) = R „. 

r —>-0 


Thus, add equations of (H21) (Il3l) and flU} together, we can give lim 

r—t 0 


dA mri (r,R mk ,e m k) 

dr 


as equation m- 


7.3 Derivation of H™~ n 

According equation (HD, H™ n can be further expressed as equation (f45]l . 


+00277 


H™~ n =\ n vf f li rn^f(e mk )f(R mk )dO mk dR mk 


+0027T 


0 0 


r —>-0 


—A n t? J J lim d ™ n \ m R rn kexp( 7r\ m R rn j^d6 rn ] £ dR rn k 


o o 


r —>-0 


For simplicity, we define h™ 1 n , h™ 2 n and /i™ 3 n in equation of (l46ll . 

^r^ ( 2 ) [- 2 ^cos(M 

h ™2 n=/(1) l _2 ^r" cos ( 0 mfe)arccos(^^ cos(0 mfe ))] 

Kr = / (1) [2^/(iZ{f) 2 -ffi cos(C)) 2 ] 

Thus 


Tjm—n 

H k 


=nKr+Ki n +Kj n \ 


(45) 


(46) 


(47) 


And Ef/i^ "], E[/i™ 2 n ] and E|7 i™ 3 71 ] can be derived as equations of (l48l) . (l49l) and (l50l) . respectively. 
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+oo 

nK~ n ] = 2X n X m v f 


/ cos {O m k)Rmk exp{-7T J2 (A i{R l i) 2 )}dOmkdRmk 


0 


7 r—arccos(min( 1 , )) 


N 

E 

i= 1 


N 


+°° rr>lb\2 „ 

= 2X n X m v J [27T^2-(-sin(0 mfc ))| 7r . » mfc ]i?mfc exp{- 7 r X] (Aj(i? b ) : 2 )}dR mk 

0 mfc 7r—arccos(min(l, )) J 

+ °° , 0 16n2 


+ 00 rylb\2 , --- JV 

=2A„A m u / [27T ( ^"^ -El - (min(l, -^)) 2 ]-R m fc expf- 7 r E (Xi(Rf ) 2 )}dR mk 
o m v n 2=1 


(48) 


— min(l, ) 

^[hT^ 71 ] = mk> 2XnXmV f J [~ 2 | ( g T ‘ l Za,iccos(-£ 1L -z)\- 7 = ! X=dzR m k exp{-7r E ( Xi(Ri b ) 2 )}dRmk 

o . ,, R mks V 1 - 2 i =1 




— min(l, —TTp) 




= 2X n X m v f 2^^[(-Vl - 2 2 arccos(^z))| . R J n +2 f 
o “‘“(i’-pr) o 

• /1 R-mk \ 


, y 1 , a ]A mfc exp{—7T E (Ai(Ji? , ) a )}dfl mfc 

v ( ^f )2 - z2 


= 2AnA m « f 2^[-7rEl - (min(l, %^)) 2 + 2 f" ]R mk expf^r E (A,(flf ) 2 )}dR 

0 V 71 0 V ^) 2 -* 2 i=1 


(49) 


— min(l, — t t- ) 
+oo v R lb ’ 


min(l, — 7 ^ ) 

v ’ R l J> ’ 


_ / p \ -)-oo /- N 

HK^} Z=COS ± mk) 2X n X m v J / [2y 1 - (-^zy-^R l *}R mk exp{—7r g(A,(i?f ) 2 )}di? mfc 




=2A„A m u / 2 ^j 2-[2 / v ^==p —exp{-7r E (Wf ) 2 )}dR mfc 

n n 2=1 


(50) 


Thus, H™ n can be derived as equation of m by adding the results of equations of C51) , (HU and (15U1) 
together. 


7.4 Proof of Corollary |T| 

According to the equation (ETT) . A™ -n = X^~ m can be proved by proving that the equation of H™~ n = 
holds. Based on the expression of H™~ n in equation (flUl) . the proof can be given by the equation 
(15T1) . where (d) follows from plugging R m k = i““( p m g m ) (e) follows from plugging z\ = Rm x ’ n z, and 

R l ^ n = ( tMj' )^x~ , (f) follows from plugging x = R n j , where R n j is the nearest distance of the nth-tier 
BSs to the origin. 
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+00 V R l rt ’ 

H?~ n = 8X n X m vJ J [ 
o o 


l-z 2 


(^) 2 - 


(%) 2 -* 


N 


+ V — 1 —]^(- R n') 2 exp{-7T £ (Ai(i?f) 2 )}di? mfc 


(d) 


1 — z 2 
Rib 

(^) 2 -Z 2 


+oo min(l, ™ ,n ) 

= 8X n X m v J J [ 
o o 

min(l.— fr —) , - 

(i) 0 , , „. + r 

SXnXmV J J [V 

0 0 v 1 

+oo min(1 ’^ ) 

Z = Z1 SXnXmV f f [ 

0 0 

{D Tjn—m 

~ H k 


+ V ( )° m expi- 71- E (M-Ri b „) 2 )}<^ 


AT 


1 — 2: 


N 




t?]^iK !n ) exp{-7r £ (E(i?“' n ) 2 )}dx 
1 2 — 1 


1-2 2 


( x \2 -2 


iV 

( *_ + V n ™i- z s - ]dz{RZ,n) 2 exp{-7T Y; (A i{R l * n ) 2 )}dx 

'Rib ' V 2=1 


(51) 


Figures 

Figure 1 - Example of donwlink HCNs with three tiers of BSs: high-power macrocell BSs (red square) 
are overlaid with successively denser and lower power picocells (red triangle) and femtocells (red circle). 
Figure 2 - The admission state and its transition of the typical UE. 

Figure 3 - The bad region when the typical UE moves to (r,0). 

Figure 4 - The average handover rates between tiers in lkm 2 region with different v. 

Figure 5 - The total handover rate (A;,) and handover rate between different tiers (A^~ 2 ) in lkm 2 region 
with different average velocity v and 2nd-tier BSs density A 2 . 

Figure 6 - The handover rates between tiers in lkm 2 region with different average velocity!; and 2nd-tier 
path loss factor a 2 . 

Figure 7 - The forward and reverse handover rates between the two tiers in lkm 2 region with different 
average velocity v, 2nd-tier BSs density A 2 and 2nd-tier path loss factor a 2 . 

Figure 8 -The CDF of residual time in the 2-tier BSs with different A 2 . 

Figure 9 -The average handover arrival rates in lfcm 2 region with different 2nd-tier bias factor (1? 2 ). 
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